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Abstract 

The goal of this paper is two-fold. The first aim is to present a detailed version of the 
computation of the two-loop renormalization of the magnetic coupling in hot QCD. The 
second is to compare with lattice simulations the string tension of a spatial Wilson loop 
using the result of our two-loop computation. 

1 Introduction 

Several breakthroughs have been made in the last few years in the computation of high 
order terms of observables in QCD at high temperature. In particular the pressure has 
been determined perturbatively up its last perturbative term namely g 6 ln(g), where g(T) 
is the QCD running coupling 1 . More than numerical results, these computations have 
given some insights for the knowledge of the convergence of perturbation theory in hot 
QCD. 

Indeed, due to infrared divergences, perturbation theory in hot QCD is limited 2 . At 
a given order in the weak coupling expansion, an infinite set of Feynman graphs will 
contribute to the thermal mean of an observable. This order depends on the observable 
or thermodynamic quantity that is computed. For example, the pressure of the gluon gas 
can be computed perturbatively up to but not including the order g 6 whereas problems 
arise at next-to-leading order for the electric Debye mass. So, despite asymptotic freedom, 
a non-perturbative part still remains in the hot QCD phase. 

Facing such a problem when one does perturbation theory and wanting to apply it in 
realistic cases renders mandatory the need for a coherent way to treat the perturbative 
expansion, taking care of infrared divergences. Whatever the approach used to get rid of 
infrared divergences, it is based on the dynamical screening of gluon field in the plasma, 
namely the Debye effect that gives a natural infrared regulator. 

We are going to work in the framework of Dimensional Reduction (D.R.) 3 ' 4 . It is based 
on the systematic construction of effective theories for each physical scales in the QCD 
plasma, namely for external momenta of order T, gT and g 2 T . This approch combined 
with 3D simulations has provided insights in phase transition in the Standard Model and 
also allows to compute the perturbative terms of the QCD pressure 1 . 

At really high temperature and large scales, within this method, the QCD plasma is 
described by a magnetostatic Lagrangian whose leading term is a three-dimensional (3D) 
Yang-Mills whose coupling is computed via perturbation theory. The important point is 
that the magnetostatic Lagrangian that describes hot QCD for external momenta of or- 
der g 2 T, contains all the non-perturbative physics in the region where the coupling is small. 



Our purpose is first to compute the effective coupling of this 3D Yang-Mills by inte- 
grating out the electrostatic Debye scale of order gT and then to use this result in order to 
investigate, thanks to lattice simulations, how our actual knowledge of the magnetostatic 
part is sufficient to encompass the whole non-perturbative high T physics in the case of 
the spatial Wilson loop. This article is composed as follows: 

Section (2) is a brief recall of the principles of D.R.. Section (3) focuses on the magneto- 
static action. Section (4) shows the steps to perform the two-loop computation. Section (5) 
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is the body of the computation with all intermediate and final result. Finally the last sec- 
tion will deal with the applications of magnetostatic physics in the QCD plasma, looking 
at an observable, the so-called Wilson Loop. 



2 Dimensional reduction perturbation theory 

Due to infrared divergences, perturbation theory in hot QCD is limited. That means 
that the analytical computation of the thermal mean of an observable is confronted with 
two main problems. The first one is to make a coherent perturbation theory until it 
breaks down. The second one is to evaluate numerically and with the most precision 
the contribution of the non perturbative sector. Indeed whatever the observable you will 
consider, the weak coupling expansion will be limited. 

For a review of effective theories and dimensional reduction see 7 . At equilibrium, a 
thermodynamic system is described by its partition function: 

Z(f3)=Tr phys exp{-f3H) (1) 

where H is the Hamiltonian of the system. In the imaginary time formalism, the latter 
can be written as a path integral: 



Z(f3) = J Dcj>Di>exp{- Jj J d z xdtC £ ) 



(2) 



with (f> a bosonic field and tp a fermionic field, with the boundary conditions 4>(x, ^r) = 
4>(x,0) and ip(x, ^) = — tp(x,0). So, technically, one is dealing with quantum field theory 
on a cylinder with a compact time direction, as in Kaluza-Klein theories. The periodicity 
of the time dimension is So, for large T, i.e. larger than the typical momentum scale 
that one considers, one expects to have an effective three dimensional theory. Due to time 
periodicity, as for massless bosonic fields, the euclidean propagator reads: 

- = 1 (3) 

p 2 p* + (27rnT) 2 ' V ' 

The massless bosons have a mass 2nnT that is going to regularize the infrared. But for 
n = 0, i.e static modes, infrared problems will still be there. Let us be more precise. 
In the case of QCD, we have: 

1 ^ 

Cqcd = -TrF^F^ + ^ ^(7m-°m + mjipi. (4) 

i=l 

As quarks have no zero Matsubara frequency, they always appear as heavy modes. 

But what is going to happen for the static modes? Part of the problem will be solved 
because for electric static mode Ao, there is Debye screening that is going to regularize the 
infrared. Indeed, like in classical plasmas, the chromoelectric field is screened in such a state 
of matter. The electrostatic potential will behave as Yukawa potential, and large distance 
contribution will be suppressed thanks to the exponential fall-off of such a potential. We 
describe below in a more quantitative way what is going to happen. 

In fact there will be different physical scales in the plasma, to which different effective 
theories will correspond. Let p be the value of the external momentum in the plasma. 

• for p ~ T, the theory is QCD for hard modes ( with n nonzero) 

• for p ~ gT, we have integrated out hard modes and obtained an effective 3D action 
for electrostatic and magnetostatic modes. 

• for p ~ g 2 T, we have integrated out electrostatic modes as they have a Debye mass 
~ gT. This leaves us with a magnetostatic Lagrangian. 
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Figure 1: Dimensional reduction at work. Integrating out each physical scale of the plasma will lead to a new 
effective theory. 



As we integrated out non-static modes, the two effectives theories are three-dimensional. 
This dimensional reduction is the same process as in Kaluza-Klein theories. The procedure 
is schematically described in Fig. 1. 
In a more precise way, we have: 

• for the hard modes: eq. (4). 

• for the soft modes, after integrating out the heavy modes, we have an SU(N) Yang- 
Mills + adjoint Higgs effective theory. Its Lagrangian reads: 

C E = Tr{D(A)A f +m 2 E TrA 2 ) + \ E {Tr(At)) 2 + 

+ \ E (Tr(A a ) 4 -^(TrA 2 ) 2 ) + ^TrF 2 +S£ E . (5) 

with as parameters g E (g,T), m E (g,T), etc ... and SCe represent higher orders 
operators. 

• for the ultra-soft modes, the effective theory is a magnetostatic Yang-Mills: 

Cm = ^TrF 2 +SC M . (6) 

with as parameter gM(gE,m E , etc..) and SCm represent higher orders operators. 

All theses Lagrangians are deduced from each other via perturbation theory. That 
means that they are known to a given order and practically one uses the truncated La- 
grangians, for example Cm — SCm, to do some computation. Of course it is important to 
know which orders are not taken into account. We will discuss this in the next section. 

The technical determination of the magnetostic Lagrangians will be discussed in the 
next chapter. We can already make some remarks: Cm is non-perturbative. The theory 
has a dimensionful coupling constant g 2 M ~ g 2 T. So if one wants to use perturbation 
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theory, the dimensionless expansion parameter is which is order 1 for \p\ ~ g 2 Tl 

Finally, once one has computed the parameters, one has to go to the lattice and use 
Cm to determine magnetic properties of the quark-gluon plasma. 
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3 Magnetostatic action in hot QCD 



As previously explained we are going to focus our attention on the magnetostatic ac- 
tion Cm- The reason is the following: the non-perturbative terms in the weak coupling 
expansion of a thermodymamic quantity or an observable come from the magnetostatic 
Lagrangian. So its study gives us some clues to study numerically and with precision 
the convergence of the perturbative series at high T. Namely, at which degree of preci- 
sion the magnetostatic Lagrangian has to be known in order to include the major part 
of non-perturbative effects? The goal of this first part is to determine this Lagrangian 
up to relative order g 2 . First we are going to recall explicitely how this magnetostatic 
Lagrangian is constructed by integrating out the Debye scale in the electrostatic action. 
Then, as the electrostatic Debye scale is integrated out via perturbation theory, we are 
going to discuss at which order we make the truncation (i.e. what is the order of 5 Cm) 
and show the computation of the magnetostatic coupling gM at two loop order. 

3.1 From the electrostatic action to the magnetostatic action 

The electrostatic action is determined by integrating out heavy modes in the full QCD 
Lagrangian, namely gluons with a non-zero Matsubara frequency and quarks. Its domain 
of validity is for external momenta of order gT. Then for momenta of order g T one has 
to integrate the Debye scale. A well-known shortcut is by making a background expansion 
of the effective potential for the magnetostatic field 8 . 

Consider Sm the full magnetostatic action containing the 3D Yang-Mills and all other 
higher order terms and Se the electrostatic action. The static fields in the electrostatic 
Lagrangian are the massless magnetostatic At field and the electrostatic Ao field with mass 
the- In order to integrate out the Debye scale, one has to introduce a background field for 
the magnetostatic field, namely, 



where the integration is from external momenta of order g 2 T and for diagrams with at 
least one Ao in the internal loops. This effective action is computed via perturbation 
theory. That means that one has to use truncation when some computations (namely 
lattice computations in this case) has to be performed with this magnetostatic action. 

3.2 Truncated magnetostatic Lagrangian and precision 

Now the aim is to see what is the relative error made by using the truncated magnetostatic 
Lagrangian. 

At this point, it is interesting to recall again that all of the non perturbative part 
is contained in the magnetostatic Lagrangian. So if one wants to know more about the 
structure of the non-perturbative terms following a perturbative expansion, for example for 
the order up to and including g for the pressure, its dominant contribution will certainly 
come from the zero-loop magnetostatic Lagrangian, namely a 3D Yang Mills theory. In fact 
the whole question is to mesure with lattice simulations the weight of this contribution and 
also to see if the high order operators in the magnetostatic Lagrangian will be important. 
This will be the aim of the second part of this paper. 

Now, let us take a look at higher order operators. One has to compute the one loop 
expansion in eq. (8). It gives two operators, namely 9 E TrF ^ F jk F kz anc j g E (DiF tJ ) 9 ^^at 



A = B + g E Q 

with B the background and Q the quantum fluctuations. 
Then formally we can construct the magnetostatic action: 



(7) 




(8) 



E 



E 
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is the order of such terms? For example, for the second term the integration of the hard 
modes have provided that at zero-loop qe ~ g 2 T and at one- loop i»e~ gT. So, as in this 
sector, D ~ g 2 T, we conclude that 8Cm is of relative order g 3 . The same occurs for the 
first term. 

To put it in a nutshell, to compute corrections of order g to the 3D magnetostatic Yang- 
Mills, the one-loop renormalization of gu is needed. For the correction of order g 2 to the 
3D magnetostatic Yang-Mills, the two loop renormalization of gu is needed. Finally for 
correction of order g 3 to the 3D magnetostatic Yang-Mills, the three-loop renormalization 
of jm but also the first term in 8 Cm are needed. Now let us go into more details for what 
is needed for the two- loop renormalization of gu- 



3.3 The coupling of the magnetostatic Lagrangian 

As we have seen, to determine 5Cm at two loop order is to determine its coupling <?m at 
two loop. 

For this we need to compute the two-point function for the background field with at 
least one Ao into the loop. We calculate the fluctuations around the background in a path 
integral: 

exp 1 -S M {B)= /DA DQ l exp4-(-S B -irr(AQ l ) 2 ). (9) 

9m J g E £ 

We added a general background gauge term. The resulting action Sm{B) is gauge invariant 
to all loop orders and the renormalization of the coupling is identified from the background 
field two point function at a momentum p = 0(g 2 T): 

exp ( - -Lsm(-B)) = exp ( - \s M {B)) (l + (ff + ff + ...)S M {B)). (10) 

9m 9e 

Here Fi is the sum of all Feynman diagrams for the two point function of the background 
field with i loops and reads: 

Fi = Fi r (S lm p 2 - pip™,)- (11) 
This leaves us with the relation, using eq.(10): 

4- = ^r-F?-F?, (12) 
9m 9e 

with 



glF? = -_L ^ (1 3 ) 

48nmE 

g\ F \r = (14) 

y 4608 V 7rm B ^ ; 
The next section shows the details of the computation that leads to this result. 



4 General and systematic method to compute diagramms 

The following section aims to recall the systematic methods in order to compute a high 
order loop expansion. It deals of course with well known facts but it is important to really 
go into detail in each of these steps in order to see how it efficiently deals with the problem, 
which is too big to be treated by hand, but simple enough to present intermediate steps 
explicitely. 

First one has to fix the gauge for the magnetostatic field in the Lagrangian. We are 
going to take the background gauge that insure the gauge independence of our result. 
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After gauge fixing, our whole electrostatic Lagrangian reads: 



C E = 



+ 



+ 



+ 



Tr(D(B + g E Q)A f + m 2 E TrAl + \ E (Tr(A 2 Q )) 
X E (Tr(A ) 4 -l(TrA 2 ) 2 ) 

^TrFij(B + g E Q) 2 + j(Di(B)Qi) 2 
8C E . 



(16) 



(15) 



(17) 



Let us recall that the vertices and propagators involving a magnetostatic field are going 
to be gauge dependent. Here are the Feynman rules with: 

• lines for Aq 

• wavy lines for A® 

• a dot for the background field 
They reads: 



p,a,i nT6~6~6~6~IT^ -p,b,j 



p,a 



p,b 





q,b 





r,c 



igfabc(Sik(pj — Tj) + Sj k (ri — qt) + 5ij(q k — Pk)) 
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9 {fadefebc ~\~ facefebd)Sij 

-ig.fabc(ri - qi) 

-igfabc($ik(Pj -rj-^) 
+3j k (ri - qt) 
+8ij{q k -p k + ^f)) 



— g 2 (fabefecd(Sik3jl — SilSjk) 
+ fadefebc(SijSkl ~ SikSjl) 
+ facefebd{SijSkl — SuSjk)) 



9 {fadefebc ~t~ facefebd)Sij 



— g 2 (fabefecd(5ik3jl — SuSjk + Sij 

+ fade febc(Sij Ski ~ SikSjl — Su-jr 
+ facefebd(Sij * — SuSjk)) 

g {fadefebc ~\~ facefebd)Sij 



The two loop graphs to be computed are shown in Fig. 2. 

The most naive way to perform the two-loop computation is by following the below 
vade-mecum. 
First: 

• contract the indices 

• scalarise integrals with no irreducible numerators using FORM 10 . 
Then: 

• Scalarise those which have irreducible numerators (to be defined later) 

The whole result can be expressed in terms of 8 scalar integrals . We have to compute 
them up to e and expand in terms of p — The result has the form: 

Ily = n(p 2 , £, d, m 2 )(p 2 S l3 - p iPj ) + Cp iPj (18) 
where following our previous convention, 

lim p ^ () n(p 2 ,(,d,m 2 ) = F? (19) 

We have to check: 

• the transversality (i.e. C = 0) 

• gauge independence 

• U.V. convergence 

• I.R. convergence 

At this point graphs 8,11,12,13,14 are zero in dimensional regularisation due to the 
presence of a massless loop. Diagrams 15 and 16 are one-loop diagrams due to the renor- 
malization of the gauge parameter 8 . 

The main body of the problem is to express the whole two-loop graphs in term of scalar 
integrals tabulated in the literature. 

5 Scalarisation, transversality and the two loop renormalisation 

5. 1 Scalarization 

The key point of the computation is to go from a set of integrals with momenta in the 
numerator to a set of "master integrals" of a few scalar integrals (8 here) to compute. We 
are going here to review the two methods of scalarisation that had been used in our two- 
loop computation. These two methods are then implemented in the Mathematica routine 
TARCER 11 . They consist of kinematical scalarisation on one hand and in a identity based 
on integration by parts on the other hand. 

Indeed, as our vertices are vectorial, after their contraction with propagators, we obtain 
for each graph the integral over internal momenta of a product of propagators times a 
numerator. The goal is then to cancel all the terms in the numerator with propagators. 
At first sight, a natural way to do this is the so-called "kinematical scalarization". It 
consists of writing the scalar products of momenta in the numerator as a combination of 
propagators. 
Let us take an example. 

Let q and r be the internal momenta, 2q.r can be written as — ((q — r) 2 +rn 2 ) + (r 2 + m 2 ) + 
(q 2 + m 2 ) — m 2 . We see here that two situations appear. If all the terms in the rewriting 
of numerators cancel with the propagators, a scalar integral will remain. This is what we 
call a reducible integral. If it is not the case, i.e. if a term in the numerator is not present 
in the denominator, the method fails. This is what we call a irreducible integral. For this 
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Graph 1 Graph 2 Graph 3 Graph 4 




Graph 5 Graph 6 Graph 7 Graph 8 




Graph 9 Graph 10 Graph 11 Graph 12 




Graph 13 Graph 14 Graph 15 Graph 16 



Figure 2: Two-loop Feynman graphs for the background two-point function. The boxes in 15 and 16 represent 
the one-loop renormalisation of the gauge parameter. 
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latter, an other algorithm must be used . Its non trivial part relies on the following 
identity 13 12 : 

n/ d% ik{ R ^ m > {q}) f\ p i-)=° (2o) 

where R is a polynomial function of internal and external momenta, L the number of loops 
in the diagram and P are the various propagators power Vj in the loop. 

Those rules have been implemented in TARCER 11 . We first use it in order to demon- 
strate the transversality of our result. The integral to appear are ( with d — 3 — 2e): 

f dd i.. 1 

J T2^F q 2 +m% 



f d d q 1 1 

J (2tt)3 ,2 +m ^ ( p+9 )2 +m ^ 






= 11 



d a q d d r 1 1 1 

(2^p'(2^F,2 + m a {p+r) 2 +m 2 {q _ r) 2 




d a q d d r 11 1 

(27r) a (2ir) 3 9 2 + m^, r' 1 +m' A E (q-r)* 



d a q d a r 1 1 1 

(27p (2ir) 3 (,!+m 2 )2 r 2 +m 2 ( 9 _ r )2 




ri a q d d r 1 1 1 

(2TF (2^F ((?+p) 2 +m; t) 2 r 2 +m2 (<j-r)2 




d"q d d r 1 1 11 1 

(2^)^ (2^F^ r2+ m 2 (,-r)'+m J (<7+p)2 (r+p)2 +m 2 



Their values are given in Appendix 1. 



5.2 Transversality for any dimension d 

The goal of this section is to check that our result is transverse. This is due to the fact 
that we compute a self energy using the background method. Before presenting the result, 
it must be noticed that only the scalarisation process is necessary to prove this. Indeed 
, one need not to express the master integrals in terms of external momenta, even not to 
specify the dimension of space. So our result is general. Let us write 

C = pillijPj = p 4 C (21) 

and check that C is zero. 
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After scalarisation, all the non transverse component can be expressed in terms of three 
master integrals. These master integrals are tabulated in Appendix 1. As defined in the 
last subsection, the propagators and vertices are all scalars . Thin lines represent massless 
propagators and bold line represent massive propagators of mass tue- For each diagram, 
the result is given in Appendix 2 in terms of the master integrals. Here is the result for 
the first graph. 

Cx = 

d - d 2 + 6£ - 5d£ + d 2 j 
12(-3 + d) 

-24m 2 + 34dm 2 - 10d 2 m 2 - 2p 2 + Sdp 2 
12(-2 + d) 

-d 2 p 2 60m 2 £ - 50dm 2 £ + 10d 2 m 2 £ 
12(-2 + d) 

10p 2 g - 7dp 2 £ + d 2 p 2 £ 
12(-2 + d) 

m 2 (4m 2 +p 2 )(l - d - 2£ + d£) 
3(-2 + d) 

(22) 

What we have learnt is that, whereas the whole result is transverse, individuals graphs 
are not. Moreover, the process of scalarization is sufficient to show explicitely transver- 
sality. It will not be the case to show gauge independance of the coupling, where the 
expressions of master integrals are needed. 

When all non transverse contributions are summed up, then the whole result is (see 
Appendix 2). The latter of course is true for any dimension. 

5.3 Scalarisation of the transverse part 

Now that we have demonstrated that our result is transverse, we want to compute the 
renormalisation of the coupling. We have to scalarize the transverse part. The complexity 
of the following computation is due to the presence of a general gauge parameter £ in the 
Lagrangian, that is reflected by the gauge dependence of the magnetic gluon propagator. 
But also, due to the presence of a magnetic background, to the gauge dependence of the 
vertices. Let us call F12 the value of a two-loop graph contracted with <5y. We have 

^Fi 2 ^p 2 Fl r (23) 

i 

For the first graph, the result reads: 



Fl 2 = 
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48m 4 - 96dm 4 + 60d 2 m 4 
12(-4 + d)(-3 + d)m 2 p 2 (4m 2 + p 2 ) 

- 12d 3 m 4 + 36m 2 p 2 - 50dm 2 p 2 
12(-4 + d)(-3 + d)m 2 p 2 (4m 2 + p 2 ) '\J 

17d 2 m 2 p 2 - 3d 3 m 2 p 2 + 12p 4 - 6dp 4 /^"N 
12(-4 + d)(-3 + d)m 2 p 2 (4m 2 + p 2 ) ( M 

-144m 4 £ + 192dm 4 g-84d 2 m 4 g ( f^\ 
12(-4 + d)(-3 + d)m 2 p 2 (4m 2 + p 2 ) V^J 

12d 3 m 4 £ - 36m 2 p 2 £ + 48dm 2 p 2 £ 
12(-4 + d)(-3 + d)mV(4m 2 + p 2 ) / 

-21d 2 m 2 p 2 C + 3d 3 m 2 p 2 £ 
12(-4 + d)(-3 + d)mV(4m 2 + p 2 ) 



4m 2 - 2p 2 + dp 2 
2(-4 + d)(-3 + d)m 2 




2m 2 +p 2 a 

2(-4+d) 

384m 4 - 720dm 4 + 408d 2 m 4 - 72d 3 m 4 
12(-4 + d) (-2 + d)p 2 (4m 2 +p 2 ) 



144m 2 p 2 - 200dmV + 134d 2 m 2 p 2 
12(-4 + d) (-2 + d)p 2 (4m 2 +p 2 ) 




-30d 3 m 2 p 2 - 32p 4 + 18dp 4 + 5d 2 p 4 - 3d 3 p 4 
12(-4 + d) (-2 + d)p 2 (4m 2 +p 2 ) 

-1152m 4 £ + 1392dm 4 £ - 552d 2 m 4 £ + 72d 3 m 4 £ 
12(-4 + d)(-2 + d)p 2 (4m 2 +p 2 ) 

-480m 2 p 2 £ + 564dm 2 p 2 £ - 222d 2 m 2 p 2 £ 
12(-4 + d)(-2 + d)p 2 (4m 2 + p 2 ) 
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+ 



30d 3 m 2 p 2 £ - 48p 4 g + 54dp 4 g - 21d 2 p 4 £ 
12(-4 + d)(-2 + d)p 2 (4m 2 + p 2 ) 

3d 3 p 4 £ _ 

12(-4 + d)(-2 + d)p 2 (4m 2 + p 2 ) 

36m 4 - 66dm 4 + 36d 2 m 4 - 6d 3 m 4 - lldm 2 p 



-e 
e- 



3(-4 + d)(-3 + d)(-2 + d)p 2 
14d 2 m 2 p 2 - 3d 3 m 2 p 2 - 6p 4 + 3dp 4 



3(-4 + d)(-3 + d)(-2 + d)p 2 

-108m 4 £ + 126dm 4 £ - 48d 2 m 4 £ + 6d 3 m 4 £ 
3(-4 + d) (-3 + d) (-2 + d)p 2 

-36m 2 p 2 £ + 48dm 2 p 2 £ - 21d 2 m 2 p 2 £ + 3d 3 m 2 p 2 £ 
3(-4 + d)(-3 + d)(-2 + d)p 2 

(-l + d)p 2 (2m 2 + 3p ; 



12(-1 + d)m 2 p 2 (4m 2 + p 2 ) 
(-l + d)p 2 




( ) 



2 






4(-l + d)m 2 

(-l + d)p 2 
2(-l + d)(4m 2 +p 2 ) 

(-1 + d)p 2 (-12m 4 + 4dm 4 - 14m 2 p 2 + 4dm 2 p 2 
12(-2 + d)(-l + d)m 2 p 2 (4m 2 +p 2 ) 

(-l+d)p 2 (2m 2 + 3p 2 ) 




6(-2 + d)(-l + d)p : 
(-l+d)p 2 







2(-l + . 




(24) 



Considering the example above, further general remarks can be made. First, compared 
with the non-transverse part of the first graph that has been shown previously, the trans- 
verse part involves more master integrals and is a little bit more complicated. Moreover, 
in the scalarisation process, some poles in e appear. They will of course cancel, witness- 
ing U.V. convergence. But attention has to be paid when one has to compute the final 
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Graph no 


P 2 Fr 


1 

2 
3 
4 
5 
6 

7 

8 

9 
10 
11 
12 
13 
14 
15 
16 


p 2 (29+24£) l+£+e(2+£-4(l+£)log(2))+4e(l+£)log(£) 
9216-n- 2 1 128£7r^ 
p 2 (20+6£) 3(2+£-2e(-l+£(-l+log(4))+log(16))+4e(2+£) log(£))) 
1536ir 2 12H€tt' a 

p 2 (-43+12£) 4+ £ (l+£-161og(2)) + 16elog(-£) 
4608-n-^ 1 64e7r 2 
_2_£ p(15+0) , P 2 (13+2f) 3(l+£) 2 +4 e (-8+£)£log(^)+4 e (3+2£(7+£))log(,£ r ) 
32p7r 1B36tt 1 3072-/T 2 768e7r 2 £ 

p(15+£(4+£)) p 2 (23+12?) 2 e +£+4 e £ log(^j) 
1536tt 4608tt 2 1 _ 128e7r 2 
o , 22-24(l+log(=£-)) 
(64 e7 r 2 ) 1 128jr 2 
6(2+«)+e(12(l+5)-p 2 (2+e))+24e(2+5)log(2£ r ) 
512e7r 2 



l+2£ | -2«-4(l+20-4(-l-2£)(l+log(^L)) 

256e7r 2 £ 1 256tt 2 £ 
16+8£ 
256p7r 






r>(l+£)(4+£) 

15367T 
p(2+£) 
768ir 


sum 


-19p 2 
4608tt 2 



Figure 3: Final result for each graph up to 0(p 2 ). 



result. Indeed some final parts when e — + are provided by the product of a i coming 
from scalarization and the first term in e of a master integral. Indeed not all the master 
integrals have to be computed up and includind 0(e). 

5.4 Result 

Once one has computed the master integrals and develop them for p 2 <C m%, the table in 
Fig. 3 is obtained. So finally 

i 

This result has of course no pole in e, is £ independant, has no pole in p, is n independant 
whereas each graph does not have these properties. This is a really strong check of our 
computation. Going into more detail will show that this check is robust in the sense that 
an error in a integral will show up through one of this checks. So the result is simply: 

Our main result, eq.(14), shows that the smallness of the corrections to the magnetic 
coupling does persist in two-loop order. In fact, at 2T C the coupling for 3 colours equals 14 
g% = 2.7 and the two-loop correction is about a third of the one loop correction (itself 
about 3 percent). 
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Our result is of importance in analyzing the purely magnetic quantities, such as the 
spatial Wilson loop, and the magnetic mass. In particular it is crucial in connecting the 
lattice results from the magnetic action to those obtained from the electric action, and 
ultimately to those of four dimensional simulations. 



6 Spatial Wilson Loop 

The goal of this part is to test the applicability of 3D physics at intermediate T. 

As has been recalled previously, the non-perturbative part in the weak coupling ex- 
pansion of an observable comes from its magnetostatic thermal average. The dominant 
term of this non-perturbative part comes from the 3D dimensional Yang-Mills with its 
coupling computed a zero-loop order. Its next term is a 3D dimensional Yang-Mills with 
its coupling computed a one-loop order. Then a three-dimensional Yang-Mills with its cou- 
pling computed a two-loop order. Then the next term is the mix of a three-dimensional 
Yang-Mills with its coupling computed at three-loop order and higher order operators. 

Our actual knowledge of the non-perturbative part is now the 3D Yang-Mills action 
with its coupling computed at two-loop order. The natural question is: does this precision 
encompass almost the whole non perturbative contribution to the thermal average of an 
observable ? 

Let us consider an observable which has its dominant contribution from the 3D Yang- 
Mills sector, unlike the pressure whose dominant contribution is the Stefan-Boltzmann 
term due to hard modes. That means that at infinite temperature, its thermal average 
is computed with the zero-loop magnetostatic Lagrangian, namely the three-dimensional 
Yang-Mills action. 

Such an observable is the spatial Wilson loop in the fundamental representation: 

W(L) = TrPexp (i <j> gA- dl). (27) 

As L is purely spatial, it measures the magnetic flux in the plasma. The thermal average 
of this spatial Wilson loop shows area behaviour with a surface tension er(T). As this is a 
purely magnetic quantity, we expect from dimensional arguments that *Ja — cg 2 M where 
c is a nonperturbative proportionality constant. Indeed, as the average of the loop is due 
to long distance correlations, hard modes will not have any effect on the thermal average. 
In the same way, for soft modes, we can integrate out the Aq field, which is what we have 
done in the previous section while constructing Lm, as the loop does not depend on Ao ! 
Finally, we have: 

J DAexp-S M (A) 

and it gives, as 5Lm is of relative order g 3 , 

a(T)=c 2 g 4 M (l + 0(g s )). (29) 

Now the aim is to fit with our formula for gu the proportionality constant c. For that, 
we have to go to the lattice. We are going to fit as a function of Indeed, gu is a 
function of gs and tue which are functions of g and T. So, for N = 3: 

2 



^ - 2 r- 1 (l + ,f- + ^) (30) 



with 



g- 2 = il(Lo S ( J) + Log{-^-) + 1.90835...). (31) 
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0.85 



Figure 4: Plot of -S= versus and fit with our two-loop formula. Each of the 10 points is shown with its error 
bar. 



The n dependance will be discussed in next section. So thanks to these formulas, one can 
fit -j- versus ^- in order to determin c and A Tc . Data points for the SU(3) spatial string 

tension has been taken from 15,16 for T > 2T C . So we have 10 points. The fit is shown in 
fig 4. 

The fit gives: 



T= = L78 ( 12 ) ( 32 ) 
c = 0.505(11). (33) 

In order to estimate the error in our fit, we have plotted for those 10 points the confidence 
levels taken from the x 2 distribution. Compared with the results in the literature our 
values are incompatible with them. 

Then we have also fitted data with T > 2.5T C . This leaves us with 6 points. We 
obtained 

= 1.57(20) (34) 



c = 0.488(18) (35) 

We have plotted for the 4 more precise of those 6 points the 95 percent confidence level 
contour taken from the \ 2 distribution in Fig. 5. Let us draw our conclusions. 

Despite the fact that the region of confidence is larger, the value of , T< = and c is still 

MS 

incompatible with refs. 17 > 18 > 19 . To see compatibility one obviously needs string data for 
higher T. It might also be recommendable to use other Monte Carlo updates, like the one 
described in ref. 20 , which is a version of the Liischer-Weisz algorithm 21 . 
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1.6 



1 . 4 



1.2 




A 



0.375 0.4 0.4250.450.475 0.5 0.5250.55 



Figure 5: Plot of various confidence levels in the plane c- for the 4 data points with the smallest error 

MS 

bar. The widespread ellipse on the left delimits the 95 percent confidence level for the parameters. The central 
values are c = 0.481 and 



1.27. The ellipse on the r.h.s. represents the actual values for these parameters 



from the literature 23,17,18 



6.1 Renormalisation scale dependence 

In the previous section, we have taken implicitely, as renormalisation scale, the one that 
cancels the one loop renormalisation due to hard modes of the electric coupling, namely 
fir = fJ-T — 47rTe- 7E eTr 22 . 

This n dependence appears explicitely when one tries to express the magnetic coupling 
qm as a function of the QCD coupling g. In fact the renormalization scale will appear in 
two ways. The first is through the renormalization of the QCD coupling, and the second 
is in the expression of qe as a function of g. 

Namely, for a pure gauge theory, the electrostatic coupling reads 

g %=g'(j i )T(l + ^-(Log{^-)) (36) 

24 /IT 

where 



^ = iJwA^) < 37 > 

If one combines those two equations, evaluating in /i = fir + <5/i, one finds that, in our case 
with no fermions, the first term in Sfi will disappear when one expand g% near fi = fir- 
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6.2 Conclusion 



We have seen that we have not yet captured the full 4D physics with our truncation. 
That means that higher order terms are needed to make the link. This has been done 
recently by considering the two loop renormalization of qb 24 • The conclusion, also with 
the fact that we need more lattice data, is that at the temperatures that we have considered, 
perturbative and non-perturbative effects are both important. The next step is to compute 
the magnetostatic Lagrangian at higher order, determining for example the leading higher 
dimensional operators. This will be the purpose of a forthcoming paper 9 
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Appendix 1: Scalar Integrals 



In this appendix, we compute the master integrals that appear in the non-transverse part. 
The computation is carried out for d = 3 — 2e. For e = 0, these integrals have been tab- 
ulated in Ref . 25 . For some of them, when a - appears in the process of scalarisation for 
one of the Feynman graphs, we have to compute the first term in the e expansion. This 
has been done using the same techniques as in 25 . At the end, one has not to forget that 
we are only interested in the result for p 2 <C m% where p is the external momentum. We 
have used the same notations as in Ref. 25 , namely, the measure of integration is J 
and thin lines represent massless propagators and bold line represent massive propagators 
of mass m_E. 



o 




-^(l + 2 £ )(l + log^-))+0(e 2 



ArcTan-rr^ — .-, n 

I -^-(2e)R 1 (p 2 ) + 0(e 2 ) 





(4*-) 



* (J-+3/2-2m E 



L °a-TTT 



-) + 0(e) 



e 





(± + 1/2 + Log^-) + 0(e) 



(It) 



ArcTan^ 

(4-7r) 2 2m eP 



2mE 1 + (-2e)R 2 (p 2 ) + 0{e 2 ) 



For the two following master integrals, only the part in e = is required. Moreover, 
they have been developed in p = 





~ 2567rp 2 m| 2567r 2 m|, + ) 

Then, the computation contained in the main body of this paper (see Appendix 3) 
shows that we have to determine the first term in e contained in two scalar integrals. 
Namely: 
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Rl(p ^ = zl L ^^ +p ^ L09 ff^- 2 + tf) (38) 
vr ; 2 87rm B ^ 96m|7r v ^ ; v ; 

1 + Log ,„ ^ , 2 -17- 12Loot^-^ 
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Appendix 2 

In this Appendix, the non-transverse part, contracted with piPj is computed for each graph. 
Ci = 




( ) 2 + 



d - d 2 + 6£ - 5dg + d 2 g 
12(-3 + d) 

-24m 2 + 34dm 2 - 10d 2 m 2 - 2p 2 + 3dp 2 
12(-2 + d) 



-d 2 p 2 60m 2 £ - 50dm 2 £ + 10d 2 m 2 £ 
12(-2 + d) 



10p 2 £ - 7dp 2 £ + d 2 p 2 £ 




e 



12(-2 + d) 

m 2 (4m 2 + p 2 )(l - d - 2£ + d£) 
3(-2 + d) 





(40) 



C 2 




( ) 2 



(-l + d)(-l + £) 



-3m 2 + 4dm 2 - d 2 m 2 - 2p 2 + 3dp 2 - d 2 p 2 + 3m 2 £ 



2(-2 + d) 



-4dm 2 g + d 2 m 2 £ + 8p 2 £ - 6dp 2 £ + d 2 p 2 ( > 



e 



+ 



2(-2 + d) 
(-l + d)m 2 (4m 2 +p 2 )(-l + £) 



e 







2(-2 + d) 

(41) 



C 3 = 

-6m 2 + 10dm 2 - 4d 2 m 2 - 6p 2 + 9dp 2 ■ 
+ 12(-3 + d)m 2 ^ " 




2 



+ -3d 2 p 2 + 18m 2 g - 18dm 2 g + 4d 2 m 2 £ ^ ^ 2 



ft , . 
12(-3 + d)m 2 
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(-18m 2 + 8dm 2 + 2p 2 - dp 2 )(l - d - 3£ + d£) 
6(-2 + d) 



m 2 (4m 2 + p 2 )(l - d - 3g + dQ 
3(-2 + d) 



9 2 + 9dp 2 - 3d 2 p 2 + 15dm 2 £ - 33d 2 m 2 £ 
24(-5 + d)(-3 + d)dm 2 £ " 



(( )) 2 




24(-5 + d)(-3 + d)dm 2 £ 

12d 2 p 2 e - 3dVg - 75dm 2 g 2 /"~"V 2 
24(-5 + d)(-3 + d)dm 2 £ \ )' 

85d 2 m 2 C 2 - 29d 3 m 2 e 2 + 3d 4 m 2 £ 2 
24(-5 + d)(-3 + d)dm 2 £) 



-21m 2 + 30dm 2 - 9d 2 m 2 + 57m 2 £ - 46dm 2 g 
12(-2 + d) 

9d 2 m 2 C + 2p 2 i - dp 2 £ 





12(-2 + d) 



m 2 (4m 2 + p 2 ) (3 - 3d - 7g + 3d£) 
12(-2 + d) 



-3 + 6d - 3d 2 + 15g - 14dg + 3d 2 £ 
12(-3 + d) 




-21m 2 + 30dm 2 - 9d 2 m 2 + 57m 2 g 
6(-2 + d) 

-46dm 2 £ + 9d 2 m 2 £ + 2p 2 £ - dp 2 £ 
6(-2 + d) 



m 2 (4m 2 +p 2 )(3 - 3d - 7£ + 3d£) 
6(-2 + d) 



C 6 = 



CV 



C 9 



Cio = 
Ci 5 = 
Cie = 



(-2 + d)(-l+d)p 2 

4(-3 + d)m 2 U ' ; 





-3m 2 + 4dm 2 - d 2 m 2 - 2p 2 
4(-2 + d) 

3dp 2 - d 2 p 2 + 3m 2 £ - 4dm 2 £ + d 2 m 2 £ 
4(-2 + d) 

8p 2 £ - 6dp 2 £ + d 2 p 2 £ 




4(-2 + d) 

(-l + d)m 2 (4m 2 +p 2 )(-l + Q 
(4(-2 + d) 



e- 



(-2 + d)(-l + d)p 2 (l-£ + dg) 
8(-5 + d)(-3 + d)dm 2 f 




( ) 2 



When all the contributions are summed up, the result is zero. 
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Appendix 3 

In this section the expression of the transverse part of each graph is given in terms of the 
8 scalar integrals: 

Fl 3 = 



48m 4 - 96dm 4 + 60d 2 m 4 
12(-4 + d)(-3 + d)mV(4m 2 + p 2 ) 

-12d 3 m 4 + 36m 2 p 2 - 50dm 2 p 2 
12(-4 + d)(-3 + d)mV(4m 2 +p 2 ) 

17d 2 m 2 p 2 - 3d 3 m 2 p 2 + 12p 4 - 6dp 




( ) 2 




( ) 



12(-4 + d)(-3 + d)m 2 p 2 (4m 2 + p 2 ) 

-144m 4 £ + 192dm 4 £ - 84d 2 m 4 £ 
12(-4 + d)(-3 + d)m 2 p 2 (4m 2 + p 2 ) 



12(-4 + d)(-3 + d)m 2 p 2 (4m 2 + p 2 ) 



-2ld 2 m 2 p 2 £ + 3d 3 mVg 
12(-4 + d)(-3 + d)m 2 p 2 (4m 2 + p 2 ) 



4m 2 - 2p 2 + dp 2 





2 




k( ) 




( ) 



_> 




(( If 




■( ) 




2(-4 + d) (-3 + d)m 2 

2m 2 + p 2 
2(-4 + d) 1 

384m 4 - 720dm 4 + 408d 2 m 4 - 72d 3 m 4 
12(-4 + d)(-2 + d)p 2 (4m 2 + p 2 ) 

144m 2 p 2 - 200dm 2 p 2 + 134d 2 m 2 p 2 
12(-4 + d)(-2 + d)p 2 (4m 2 + p 2 ) 

-30d 3 m V - 32p 4 + 18dp 4 + 5d 2 p 4 - 3d 3 p 4 
12(-4 + d)(-2 + d)p 2 (4m 2 + p 2 ) 

-1152m 4 C + 1392dm 4 £ - 552d 2 ra 4 £ + 72d 3 m 4 £ 
12(-4 + d)(-2 + d)p 2 (4m 2 + p 2 ) 
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-e 
e- 





-480m 2 p 2 £ + 564dm 2 p 2 £ - 222d 2 m 2 p 2 £ 
12(-4 + d) (-2 + d)p 2 (4m 2 +p 2 ) 

30d 3 m 2 p 2 g - 48p 4 £ + 54dp 4 g - 21d 2 p 4 £ 
12(-4 + d) (-2 + d)p 2 (4m 2 +p 2 ) 

3rfVe _ 

12(-4 + d)(-2 + d)p 2 (4m 2 +p 2 ) 
36m 4 - 66dm 4 + 36d 2 m 4 - 6d 3 m 4 - lldm 2 p 



e 





3(-4 + d)(-3 + d)(-2 + d)p 2 
14d 2 mV - 3d s m 2 p 2 - 6p 4 + 3dp 4 



3(-4 + d)(-3 + d)(-2 + d)p : 

-108m 4 £ + 126dm 4 £-48d 2 m 4 £ + 6d 3 m 4 £ _ 
3(-4 + d)(-3 + d)(-2 + d)p 2 

-36m 2 p 2 £ + 48dm 2 p 2 £ - 21d 2 m 2 p 2 £ + 3d 3 m 2 p 2 £ 
3(-4 + d) (-3 + d) (-2 + d)p 2 

(-l + d)p 2 (2m 2 +3p 2 ) y \ a 



12(-1 + d)m 2 p 2 (4m 2 + p 2 ) 
(-l + d)p 2 




( ) 





e 



4(-l + d)m 2 

(-l + d)p 2 
2(-l + d)(4m 2 +p 2 ) 

(-1 + d)p 2 (-12m 4 + 4dm 4 - 14m 2 p 2 + 4dm 2 p 2 ) 
12(-2 + d)(-l + d)m 2 p 2 (4m 2 +p 2 ) 

(-1 +d)p 2 (2m 2 +3p 2 ) 
6(-2 + d)(-l + d)p 2 



(-l+d) P 2 
2(-l + d) 
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F2 2 = 



-18m 2 + 14dm 2 - 6p 2 + 3dp 2 ( /^~\ 
8(-3 + d)m 2 (4m 2 +p 2 ) \ J 



3(-4m 2 + 4dm 2 - 2p 2 + dp 2 ) 
8(-3 + d)m 2 




-18m 2 + 14dm 2 - 10p 2 + 5dp 2 - 48m 2 £ + 24dm 2 £ - 12p 2 £ + 6dp 2 £ 
4(-2 + d)(4m 2 +p 2 ) 




-18m 2 + 14dm 2 - 6p 2 + 3dp 2 
4(-3 + d)(-2 + d) 




3(-l + d)p 2 (2m 2 +3p 2 ) 
12(-1 + d)m 2 p 2 (4m 2 +p 2 )2 




3(-l + d)p 2 
4(-l +d)2m 2 



o 




3(-l + d)p 2 
2(-l +d)2(4m 2 +p 2 ) 



e 



3(-l + d)p 2 (-12m 4 + 4dm 4 - 14m 2 p 2 + Adm 2 p 2 ) 
12(-2 + d)2(-l + d)m 2 p 2 (4m 2 + p 2 ) 




3(-l + d)p 2 (2m 2 +3p 2 ) 
6(-2 + d)2(-l + d)p 2 




3(-l + d)p 2 
2(-l + d)2 




F3 2 = 



+ 



(-2 + d)(-4m 2 + 6dm 2 - 2d 2 m 2 - 4p 2 + 5dp 2 - d 2 p 2 + 12m 2 £ - 10dm 2 £ + 2d 2 m 2 £) 

(4(-4 + d)(-3 + d)m 2 p 2 



(-1 + rf) 2 
2(-3 + d) 





2G 



+ 



(-16m 2 + 6dm 2 + dp 2 )(l - d - 3g + dQ 
2(-4 + d)p' 2 

2m 2 (-6m 2 + 2dm 2 - 2p 2 + dp 2 )(l - d ~ 3£ + d£) 
(-4 + d)(-3 + d)p 2 



F4 2 



+ 



(-2 + d)(240m 2 p 2 - 240dm V - 20160m 4 £) 
1920(-5 + d)(-3 + d)m 4 p 2 £ 






(-2 + d)(+52832dm 4 £ - 37440d 2 m 4 £ + 9696d 3 m 4 £) 
+ 1920(-5 + d)(-3 + d)m 4 p 2 £ 



(-2 + d)(-832d 4 m 4 £ + 1864m 2 p 2 £ - 4008dm 2 p 2 £) 
+ 1920(-5 + d)(-3 + d)m 4 p 2 £ 



(-2 + d)(+3352d 2 m 2 p 2 £ - 1208d 3 m 2 p 2 £ + 128d 4 m 2 p 2 £) 
+ 1920(-5 + d)(-3 + d)m 4 p 2 C 

(-2 + d)(-420p 4 £ + 509dp 4 £ - 225d 2 p 4 £ + 43d 3 p 4 £) 
+ 1920(-5 + d) (-3 + d)m 4 p 2 £ ' 

(-2 + d)(-3d 4 p 4 g + 96320m 4 £ 2 - 115104dm 4 g 2 ; / v 
+ 1920(-5 + d)(-3 + d)m 4 p 2 C ' ' 






(-2 + d)(+48768d 2 m 4 £ 2 - 8800d 3 m 4 £ 2 + 576d 4 m 4 g 2 ) 
+ 1920(-5 + d) (-3 + d)m 4 p 2 C 



(-2 + d) (-3128m 2 p 2 ^ 2 + 5832dm V^ 2 - 3432d 2 m 2 p 2 £ 2 ) 
+ 1920(-5 + d)(-3 + d)m 4 p 2 £ 



(-2 + d)(+792d 3 m 2 p 2 C 2 - 64d 4 m 2 p 2 £ 2 - 60p 4 £ 2 
+ 1920(-5 + d)(-3 + d)m 4 p 2 £ 

(-2 + d)(+107dp 4 £ 2 - 59d 2 p 4 £ 2 + 13d 3 p 4 £ 2 - d 4 p 4 £ 2 




+ 



1920(-5 + d)(-3 + d)m 4 p 2 £ 

1680m 4 - 6840dm 4 + 6480d 2 m 4 - 1320d 3 m 4 
1920m 4 

27 




O 



-240mV + 640dmV - 540d 2 m 2 p 2 + 140d 3 m 2 p 2 
+ 1920m 4 

56p 4 - 66dp 4 + 25d 2 p 4 - 3d 3 p 4 - 16080m 4 g + 23640dm 4 £ 
+ 1920m 4 

-8400d 2 m 4 £ + 840d 3 m 4 £ + 400m 2 p 2 £ - 720dm 2 p 2 £ 
+ 1920m 4 

380d 2 m 2 p 2 £ - 60d 3 m 2 p 2 £ + 8p 4 £ - 14dp 4 £ 



o 



+ 



+ 



1920m 4 

7d 2 p 4 £ - dVC - 6480m 4 g 2 + 9720dm 4 £ 2 
1920m 4 





-3600d 2 m 4 £ 2 + 360d 3 m 4 £ 2 - 2160m 4 f 
+ 1920m 4 

3240dm 4 £ 3 - 1200d 2 m 4 £ 3 + 120d 3 m 4 £ 3 
+ 1920m 4 

16128m 4 - 37024dm 4 + 18272d 2 m 4 - 2496d 3 m 4 




o 



+ 



c 



960m 2 p 2 




-56mV + 784dm 2 p 2 - 1336d 2 m 2 p 2 + 384d 3 m 2 p 2 
960mV 




168p 4 - 191dp 4 + 72d 2 p 4 - 9d 3 p 4 - 77056m 4 £ 
960mV 




67040dm 4 £ - 18912d 2 m 4 g + 1728d 3 m 4 £ 
960mV 




-968mV£ - 1312dm 2 p 2 £ + 1112d 2 mV£ 
960m 2 p 2 




-192d 3 m 2 p 2 g + 24p 4 g - 41dp 4 £ + 20d 2 p 4 £ - 3d 3 p 4 g 
960mV 




(4m 2 +p 2 ) (-2016m 4 + 3872dm 4 - 832d 2 m 4 ) 
960mV 
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(4m 2 +p 2 )(64m 2 p 2 - 208dm V + 128d 2 mV) 



+ 



+ 







960m 2 p 2 


(4m 2 


+ P 2 )(- 


-28p 4 + 19dp 4 - 3dV + 9632m 4 £) 






960m 2 p 2 


(4m 2 


+ P 2 )(r 


-4768dm 4 £ + 576d 2 ra 4 £ - 208mV£) 



960m 2 p 2 

(4m 2 +p 2 ) (272dm - 64d 2 m 2 p 2 ^) 



960m 2 p 2 

(4m 2 +p 2 )(-4p 4 g + 5dp 4 £ - d 2 p 4 Q 
960m 2 p 2 



F5 2 = 



416m 2 - 576dm 2 + 264d 2 m 2 - 40d 3 m 2 




48(-3 + d)m 2 p 2 
+48p 2 - 53dp 2 + 18d 2 p 2 - d s p 2 - 608m 2 £ A \ 2 



48(-3 + d)m 2 p 2 

768dm 2 £ - 312d 2 m 2 £ + 40d 3 m 2 £ - 6p 2 £ 
48(-3 + d)mV 



„2c Rj2 m 2 c | j3„2t 

'-([ )T 





( ) 



2 




48(-3 + d)m 2 p 2 

(-1 + d)(4m 2 + p 2 )(-l + Qj-d - 4g + dQr[(-l + d)/2] 
32r(l/2 + d/2) 

8(16m 2 - 8dm 2 + 3p 2 - 2dp 2 - 24m 2 £) 
32 

8(8dm 2 £ - 6p 2 £ + 2dp 2 £) 
32 





+ 3(-l + d)m 2 (-l + £)(-d - 4£ + d£)r[(-l + d)/2] 

29 



o 



192m 2 - 186dm 2 + 42d 2 m 2 + 26p 2 - 15dp 2 



48m 2 



d 2 p 2 - 264m 2 g + 204dm 2 g - 36ri 2 m 2 £ 
48m 2 



-2p 2 £ + 3dp 2 £ - d 2 p 2 i - 24m 2 : £ 2 + 30dm 2 £ 2 
48m 2 



-6d 2 m 2 £ 2 
48m 2 



-4m 2 + 4dm 2 - 3p 2 + 2dp 2 
8(-3 + d)m 2 




o 2 (8m 2 + p 2 ) 




1664m 2 - 2096dm 2 + 872d 2 m 2 - 120d 3 m 2 
24(-2 + d)p 2 



288p 2 - 283dp 2 + 74d 2 p 2 - 3d 3 p 2 - 2432m 2 £ 
24(-2 + d)p 2 



2768dm 2 £ - 1016d 2 m 2 £ + 120d 3 m 2 £ - 114p 2 £ 
24(-2 + d)p 2 



121dp 2 £ - 38dV £ + 3dV £ 
24(-2 + d)p 2 

2496m 4 - 3040dm 4 + 1216d 2 m 4 
24(-3 + d) (-2 + d)p 2 




160d a m 4 + 960m 2 p 2 - 1068dm 2 p 2 + 376d 2 m 2 p 2 
24(-3 + d) (-2 + d)p 2 



'Ud 3 rn 2 p 2 + 96p 4 - 83dp 4 + 18d 2 p 4 - d 3 p 4 
24(-3 + d) (-2 + d)p 2 



-3648m 4 £ + 4000dm 4 £ - 1408d 2 m 4 £ + 160d 3 m 4 £ 
24(-3 + d) (-2 + d)p 2 
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-936mV£ + 1044dmV£ - 376d 2 mV£ 
24(-3 + d)(-2 + d)p 2 




44d 3 m 2 p 2 £ - 6p 4 C + lldp 4 £ - 6d 2 p 4 £ + d 3 p 4 £ 
24(-3 + d)(-2 + d)p 2 




F62 = 



(-2 + d)(-l + d)d i^^V a 
4(-3 + d)m 2 l \ )' 



F7 2 = 



3(-l + d + Q 
4 




F9 2 = 

(-2 + d)(-l + d)(l - g + dg) i^"^ 
8(-5 + d)(-3 + d)m 2 £ \y' 

FIO2 = 

(-1 + d)(-14 + 43d - lid 2 + 134£ - 63d£) 
16 



(-1 + d)(7d 2 £ + 54£ 2 - 27d£ 2 ) 
16 

(-1 + d)(3d 2 £ 2 + 18£ 3 - 9dg 3 + d 2 g 3 ) 
16 
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